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ABSTRACT 

Particle acceleration at non-relativistic shocks can be very efficient, leading 
to the appearance of non-linear effects due to the dynamical reaction of the 
accelerated particles on the shock structure and to the non-linear amplifica- 
tion of the magnetic field in the shock vicinity. The value of the maximum 
momentum p max in these circumstances cannot be estimated using the clas- 
sical results obtained within the framework of test particle approaches. We 
provide here the first attempt at estimating p max in the cosmic ray modified 
regime, taking into account the non-linear effects mentioned above. 

Key words: acceleration of particles - shock waves 



1 INTRODUCTION 

Since the beginning of the investigations on the mechanism of particle acceleration at shock 
fronts, estimating the maximum energy of the accelerated particles has been one of the main 
goals, in the perspective of assessing the ability of the mechanism to explain the observed 
energies in cosmic rays. This issue has played a crucial role especially in the context of 
the so-called supernova paradigm for the origin of galactic cosmic rays. In fact, it was 
understood very soon, in the case of supernova remnants' shocks, that acceleration would 
not be efficient enough and provide sufficiently high energies unless magnetic scattering 
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could be self-generated by the particles themselves through their motion in the upstream 
plasma. The mechanism thr ough which this generation of waves could occur was identified 



i n the s treaming instability (jBelll (1l978al )) and widely discussed later by lLagage fc Cesarsky 
(jl983al ibl). In these latter papers, in particular, the maximum energy was estimated for 
shocks in supernova remnants, already believed to be candidate sources of galactic cosmic 
rays. The result of this investigation was however that the maximum energy, even in the 
presence of self-generated scattering agents, is of the order of ~ 10 3 — 10 4 GeV, quite lower 
than the energy at which the knee is observed. Today we know that this energy is also lower 
tha n the maximum e nergy observed in the proton component by the KASCADE experiment 



see 



Horandell (120051 ) for a comprehensive re yiew of expe rimental results in the knee region 



that seems to be as high as ~ 10 7 GeV (see iBlasil ( 120051 ) for a recent review). 

More recently the problem has been complicated even more by the understanding that 
particle acceleration can be rather efficient and break the test particle regime that old cal- 
culations were based upon. This important res ult was obtained independe ntly within com 
pletely d i fferen t approaches, two- fluid models (IDrurv fc V olk 



(Markov 



( Jones fc Ellison (1991); 



(119971 ): iMalkov. Diamond fc Volkl (I2000h : 



Blasil (12002 



1980 



1981)). kinetic models 



2004T ): lAmato fc Blasil (12005 . 



2006) ) and numerical approaches, 



Bel 



u sing both Monte Carlo and oth e r sim u 



( 1987 1 : Ellison. Mobius fc Paschmannl (11990 ): 



200.4 ): 



ation procedures 



Ellison, Ba ring fc Jones 



1995 



1996|); 



Kang fc Jones 



(1997 



Kang. Jones fc Gieseler! (120021 ): 



Vladimirov. Ellison fc Bykov 



( 120061 )). The non-linear dynamical reaction of the accelerated particles introduces a positive 
reaction that results in a flattening of the high energy spectra, which are no longer power 
laws, thereby enhancing the role of particles at the highest energies. Moreover a gradient in 
the velocity of the upstream fluid is produced by the pressure of the accelerated particles. 
The value of p ma x is important to determine the spectrum, and the spectrum determines 
the amplification o f the mag netic field which in turn determines p ma x • The formalism of 
Lagage fc Cesarskyi (ll983a! Jbl) is not appropriate to describe this type of situation. 

In this paper we describe the first calculation of the maximum momentum for cosmic 
ray modified shocks, including all the relevant non-linear effects. The description of the 
acceleration process, with the dynamical reaction of cosmic rays and the magnet i c field 



amp 



2006|). 



ification taken into account, is based on the recent papers by 



Amato fc Blasil (12005 . 



The paper is organized as follows: in SJ2] we describe the method that allows us to com- 
pute the acceleration time as a function of the particle momentum; in we describe the 
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procedure used to determine p ma x in the case of modified shocks and discuss the underlying 
assumptions; in §4] we present the results we find for the maximum achievable momentum 
in supernova remnant shocks; finally a critical discussion of the method adopted and a 
summary of the main results is provided in § [5l 



2 ACCELERATION TIME FOR MODIFIED SHOCKS 

The full description of the calculations that allow us to compute the spectrum of the accel- 
erated particles, the dynamics of the upstream and downstream plasmas and the strength 
of the amplified magnetic field due t o cosmic rays' s t reaming ins tability in the stationary 
regime can be found in the papers by lAmato fc Blasil (120051 . 120061 ) . 

In this section we present the main achievement of our work, which consists of a deriva- 
tion of the acceleration time t acc (p) (as a function of the momentum p) for particles accel- 
erated at a shock with an arbitrary level of shock modification. The expression we derive 
i s a generalization o f the w ell known formulae found in the context of quasi-linear theory 
( Lagage fc Cesarsky ( 1983aJ lbl) ) and reduces to them when the non-linear shock modification 
is negligible. 

In general, the transport of the accelerated particles is described by the time dependent 
diffusion-convection equation: 



dt dx 



d_ 

dx 



D 



dx 



1 du df ^ , . 
+ o^-Ptt- + Q{x,p), 
6 dx op 



(1) 



where f(t, x, p) is the distribution function of the accelerated particles, D(t, x, p) is the spatial 
diffusion coefficient and Q is the injection term, taken in the form Q(x,p) = Qo{p)S(x), with 
the assumption that injection takes place only at the shock surface on the downstream side 
(although this assumption is not strictly needed). The x— axis goes from upstream infinity 
(x = — oo) to downstream infinity (x = +oo) and the shock is at x = 0. We will use the 
index 1 (2) to refer to the region immediately upstream (downstream) of the shock. 

The problem of determining the time scale associated with the change of the distribu- 
tion function around the maximum momentum is impossible to face analytically in all its 
complexity, as implied by the non-linearity of the system: in general, /, u and D all depend 



on time, and change on comparable time sea 



on all other quantities ( lAmato fc Blasil (12005 



e s. Mo reover, / depends in a complicated way 



2006|)). 



Let us assume that we know an approximation of the solution / in the form of a solution 



of the stationary problem for a given value of the maximum momentum p max , as can be 
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calculated following lAmato fc Blasil (1200 



20061 ) . This value of p ma x can be considered as a 



guess value to start with. In this case the velocity profile u(x) can be considered as determined 
in first instance from the guess function f as de r ived b y the assumption of quasi-stationarity, 
using the calculations of lAmato fc Blasil (120051 . 120061 ). Within this approximate framework, 
we can introduce the Laplace transform of the distribution function with respect to time: 



g(s,x,p) 



dt e st f(t,x,p), 



(2) 



so that the transport equation becomes: 
d 



dg 

sg + u— = — 

ox ox 



ox 



ldu dg 1 
6 dx op s 



(3) 



Integration of Eq. [3] in the vicinity of the shock (namely between x = and x = 
leads to the following boundary condition at the shock: 

1 



jjdg 




jjdg 


dx 


2 


dx 



1 dg (s,p) 
+ 3 {U2 - Ul)P ^p— 



~-Qo(p) = o, 



(4) 



where U\ = w(0 + ), go(s, p) = g(s, x = 0, p) and we have used the continuity of the distribution 
function (at any time) at the shock. 

We can find another condition on the distribution function by integrating the transport 
equation (Eq. [3]) in the upstream section from x = —oo to x = 0~, so as to obtain: 



sl(s,p) + u 1 g (s,p) 



,du(x) 
dx g{s,x,p)— — 

-oo dx 



D 



dg_ 

dx 



1 d r° 



i 3^ dp 



dx g(s, x,p) 



du(x) 



dx 



where I(s,p) = J^dx g(s,x,p). Following Ifjlasil ( 20021 ) we introduce the quantity 
1 

dx ) 



U 



p {s,p) = Ui 



9o 



°" du 
dx g(s, x,p)- 

-oo 



which allows us to rewrite Eq. [5] as 
dg 



D 



dx 



sl(s,p) +90 u p (1 + 3^ 



1 dg 



(5) 



(6) 



(7) 



Clearly in this equation u p depends on time through the functions / and u, which again 
reflects the non-linearity of the underlying equation. 

The transport equation in the downstream fluid can be considerably simplified if we 
assume that the fluid velocity is independent of location u(x) = w 2 an d that the diffusion 
coefficient is also independent of x. In this case, for any x ^ Eq. [T] becomes: 
d 2 f 



dt dx 



D-, 



dx 2 



(8) 



If we use again the Laplace transform of the distribution function with respect to time we 
have 



dg d 2 g 

Sg + U 2 T- = D 2 — : 

ox ox z 
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(9) 



whose solution has the form g(s,x,p) = go(p)e l3x (in the downstream region x > 0) with (5 
easily found by substitution, so that: 



g(s,x,p) = g (p)exp 



This implies that 



u 2 x 
2D~? 



1 



1 + 



D 



dg_ 

dx 



:#o(p)«202, 



a 2 = 1 



AsD 2 



1 



1/2' 



(10) 



AsD 2 



1/2 



(11) 



At this point, substituting Eqs. [71TTT1 into the jump condition Eq. H]we obtain the following 
equation: 



-u 2 a 2 — Up 1 + 



1 d In u r 



3 d\np 



As already proposed by 



Malkov 



+ -(w 2 - u v)V-^ ~ si + -Qq 



0. 



(12) 



(119971 ) and lAmato fc Blasil (120051 ) , it is useful to rewrite the 
transport equation in an implicit form by integrating it between x = — oo and a generic point 
x upstream. We repeat this procedure here on the Laplace-transformed transport equation 
(Eq. [3]). This leads to the following implicit expression for g: 



1 d\ng(s, x", p) \ 



x ' „ u(x 
dx 



D(x",p) 



x 



dx"g(s, x",p) 



du I 

dx 7 ' \ + 3 



dlnp 



dx"g(s, x",p) 



(13) 



As already discussed by iMalkovl (119971 ) and lAmato fc Blasil (120051 ) , one can approximate 
this complex implicit solution as 



g{s,x,p) = g (s,p)exp 



q(p) 



dx' 



u[x 



D(x',p) 



+ 



+ sexp 



dx' 



U[X 



D(x',p) 



dx' 



D(x',p) J-oo 



dx" g (x" 



(14) 



where q(p) is now the slope of the distribution function at the shock. This e xpression is ob 
tained by expanding the previous expression for small values of Jn dx'^jr. In 



Amato fc Blasi 



D(x') ■ 

( 120051 ) it was shown a posteriori that this is a good approximation of the spectral distribution 
of the particles as a function of the spatial coordinate. 

The agreement was found to be worse when the shock is quasi-linear, namely for small 
modifications. The reason for this is in the fact that Eq. [14] satisfies the jump condition in 
the stationary regime (s — > 0) only when u 2 — > 0. A better approximation for g is in fact of 
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the form, 

g(s,x,p) = g (s,p) exp <— -11 / dx — — 

[3 V uiJ Jo D[x 



+ s exp 



dx' 



u[x 



D(x',p) 



dx' 



D(x',p) J-oo 



P) j 

dx" g {x' 



+ 



(15) 



which satisfies the jump condition exactly. Note that Eq. [15] differs from Eq. [TH only in 
that at given momentum the particle distribution is slightly more extended spatially. In the 
following we adopt Eq. [15] as an approximation of the spatial distribution of the accelerated 
particles. In evaluating the term sl(s,p) in Eq. [T^] we retain only the first term of g in 
Eq. [151 since in the end we will be interested in the limit s — > (large times). This means 
that Eq. [12] becomes: 



9o 



v 2 
where 

A(p) = 



u 2 o 2 - sA(p) - u p 1 



1 d In Up 
3 a! hip 



1/ \ dg 1_ 

+ ^2-u p )p— + -Q 



0. 



dxexp 



q(p) 



1 - 



U 2 



U\J JO 



dx 



, u{x') 



D(x' 



(16) 



(17) 



The solution of this equation is easily found to be: 



/ \ 1 V n gas,l 



3R 



sub 



s 47rpf • RtotUpip) 



exp 



p dp' 3U p (p')R 



tot 



Pinj P' RtotU p (p') 



X 



v dp' 3 
Xex PW T7o" 



°2 



p' 2 R to tU p (p') - 1 
or, in a more compact form, 

9o(s, p) = -K{p) exp [h(p, s)} , 



x exp 



p dp' 3sA(p') 
in] p' RtotUpip') - 1 



where 



K(p) 



and 



V n gas,l 



3R 



sub 



^Pinj RtotU p (p) 



exp 



p dp' 3U p (p')R tot 

Pinj P' RtotUp(p') ~ 



J Pinj P 



a 2 



3R 



tot 



sA(p') 



2 R tot U p (p') - 1 u R tot U p (p')-l 
with U p (p) = u p /uq and u = u(—oo). 



(18) 



(19) 



(20) 



(21) 



At this point we follow the elegant procedure discussed by iDruryl (1l983l ). We need to 
calculate the inverse Laplace transform 



] r+ioo g* s 

f(t,x,p) = — ds — K(p) 

2,111 J-ioo s 



,h(p,s) 



We introduce the function 



0(p,t) 



1 

2vri 



+«oo 



ds e ts e h{p ' s) . 



(22) 
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One can easily show that the inverse Laplace transform of the function 

f(p,t) = K(p) f dt' <j>{p,t') (23) 
Jo 

is exactly the function go as written in Eq. [191 It follows that the solution of our problem 
is in fact / = /. In the approximate f ramework in which we are, the function <p, as in the 
linear case discussed by iDruryl (119831 ) . represents the probability that a particle starting 
with energy p in j has reached the momentum p in a time t (<p is normalized to unity). The 
mean time for the acceleration to momentum p is therefore < t(p) >= J °° dt t <p(p,t). One 



can use the following obvious property of the function 

dt <f)(p, t) e~ ts = exp [h(p, s)] , 



(24) 



to determine < t >, by deriving Eq. 
this follows that: 



< t > 



Oh 
ds 



3R 



tot 



s=0 



117 



p dp' 
i P' 



with respect to s and calculating it at s = 0. From 
RtotD 2 (p') u A(p') \ 



+ 



(25) 



RtotU p (p>) - 1 R t otU p (p') - 1 J ' 
For an unmodified shock, with spatially constant fluid velocities upstream and downstream, 
one has that q(p) = 3u±/ (ui—112) and the previous formula easily reduces to the well known: 



< t > 



El 

u 2 



£1 



(26) 

Ui - U 2 L U 2 Ui J 

if the diffusion coefficient is in the form of Bohm diffusion. For an arbitrary form of the 
diffusion coefficient the expression in Eq. is correct within a factor of order unity. 



3 CALCULATION OF P MA x 



In general the value of p ma x is determined by the competition between the acceleration time 
and the shortest among the finite age of the source, the time scale for energy losses, and the 
time scale for escape of particles with momentum p max from upstream infinity. 

Here we concentrate on the self-consistent calculation of the maximum momentum in 
those cases when it is determined by the finite age of the source (this is expected to be 
the case for supernova remnants). In such a context each val ue of p mnx there f ore reflects a 



different age r of the system. In this sense, our previous work (jAmato fc Blasil ( 12005 



2006|)) 



was based on the implicit assumption that a quasi-stationary regime could be reached on 
time scales shorter than r, but no specific recipe for the calculation of p max was adopted. 
In other work on non-linear diffusive particle acceleration at shock waves, p m ax was either 
fixed or determined according with a simple recipe based on test-particle theory. 
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Our determination of the maximum momentum proceeds through 2 different steps: 1) 
we first determine the acceleration time t acc (p) (Eq. [25]) for a given value of p%} ax (at the 
first step this is just a guess value p^J; 2) we compare t acc (p!£l x ) with the age of the 
system r and determine an updated value of p max , say p^x 1 ■ The procedure is iterated until 
tacciPmax) = r , where Pmax is now the actual maximum momentum for our problem. It is 
crucial to keep in mind that due to the non-linearity of the system, changing p max causes the 
all system to change (velocity profile, spectrum of accelerated particles, efficiency of particle 
acceleration and also the diffusion coefficient, if the scattering is due to self-generated waves). 

The results obtained by applying this procedure are illustrated in the next section. 



4 RESULTS 



Following the method outlined in § [2] we computed the acceleration time as a function of 
the particle momentum p for a system with r = 1000 yr, under different assumptions on the 
shock Mach number at infinity M and on the diffusion properties of the system. The results 
of these calculations are summarized in Fig. [TJ We assumed a velocity of the shock wave of 
5 x 10 3 km/s and considered values of Mo ranging from 5 (unmodified shocks) to 200 (already 
in the asymptotic regime for strongly modified shocks). The injection momentum was taken 
to be £ ptherm with ptherm corresponding to the maxwellian temperature T 2 and £ = 3.5. We 
also related the inj ection efficiency r\ to the c ompr ession ratio at the subsh ock through the 



recipe propo sed by iBlasi. Gabici fc Vannonil (120051 ) and already adopted by 



Amato fc Blasi 



( 120051 . 120061 ). In the top panel of Fig. [I] we plot the acceleration time up to a momentum 
p for D(x,p) = D sg , namely with self-generated scattering, with D sg = (4/37r)w£/ T the 
result of quasi-linear theory for scattering induced by self-generated Alfven waves. Here v 
and tl are the particle velocity and Larmor radius (in the background field -Bo) respectively, 
and T is the dimensio nless wave energy density per unit logarithmic band-width (see e.g. 



Amato fc Blasil (120061 ) for details). In the bottom panel of Fig. [I] we show the acceleration 



time for D(x,p) = Dsohmi with Dsohm referring to Bohm-like diffusion in the amplified 
magnetic field 5B found by integrating B\T in /c-space: Dsohm, = (i/^)vri(B /SB). In all 
of our calculations the background magnetic field was taken to be Bq = 10/iG. 

The maximum momentum obtained from our iterative procedure is seen as the intersec- 
tion with the horizontal line r = 1000 years, assumed to be the age of the system. 

Each of the curves in Fig. [1] was obtained by following the procedure outlined in § |3lThe 
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Figure 1. We plot the time (in years) it takes for the system to accelerate particles up to a momentum p (in units of mc) under 
different conditions. The different panels refer to different choices of the diffusion coefficient. In the top panel the diffusion 
coefficient is the self-generated one, while in the bottom panel we assume Bohm diffusion in the amplified magnetic field (see 
text for details). The different line types refer to different values of the shock Mach number: Mo = 5 for the solid line, Mo = 20 
for the dotted line, Mo = 100 dashed, Mq = 200 for the dot-dashed line. 



first thing to notice by comparing the two different panels in Fig. [T] is that the maximum 
achievable momentum, for a given Mach number, is roughly the same (within a factor of a 
few (2-4) at most) for the two diffusion models adopted. This observation is very important 
given the fact that the diffusion coefficient is basically unknown and not easily obtained 
from observations. 

Still the comparison between the two panels allows to highlight what the main effect of 
the different diffusion models is: in the case of strongly modifie d shocks the spe c tra b ecome 
very flat at high energies as discussed in detail elsewhere (e.g. lAmato fc Blasi ( 2006 )) and 
most of the power of the Alfven waves excited by the streaming cosmic rays is deposited 
in long wavelength waves (the ones that are resonant with high energy particles), with a 
relatively smaller fraction going into waves that are able to efficiently scatter the low energy 
particles. This is the reason for the behaviour of t acc (p) observed in the upper panel, which 
may seem surprising at first sight for the flatness of the curves referring to Mq = 100 — 200. 
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By comparing the two panels, it is clear that things change drastically, with the dependence 
of t acc on p becoming much stronger, if the energy of the waves is redistributed uniformly 
between different wavelengths, which is what one is effectively assuming by adopting a 
Bohm-like diffusion coefficient in the amplified magnetic field. 

Finally, let us discuss the behaviour of the maximum achievable momentum as a function 
of the Mach number. This is qualitatively the same for both models of diffusion: p max initially 
increases quickly with the Mach number as a result of field amplification in modified shocks 
and then decreases somewhat, although slowly, when Mq is further increased. This latter fact 
is clearly related to the more strongly modified velocity profile that leads to lower and lower 
values of the flow velocity in the vicinity of the shock and as a consequence increasingly 
longer residence times upstream. The dependence of the maximum achievable momentum 
on the shock Mach number is illustrated more clearly in Fig. [2J where we plot p ma x as a 
function of Mq for four different models of diffusion. The plot refers again to a system with 
t = 10 3 yr and all other values of the relevant parameters are the same as above, including 
the background magnetic field strength and the recipe for Pi n j. 

The solid curve refers to D(x,p) = D sg (x,p) as described above, the dashed curve is for 
D(x,p) = D Bohm (x,p), the dotted curve is for D(x,p) = D Bohm (0,p) (namely the diffusion 
coefficient equals the Bohm diffusion calculated in the amplified magnetic field at the shock 
surface, without spatial variation) and finally the dot-dashed curve refers to Bohm diffusion 
in the background magnetic field. The first aspect that clearly appears from the figure is the 
enhancement of p ma x that arises from magnetic field amplification. For high Mach number 
shocks the magnetic field strength at the shock is amplified by a factor of order 25 and 
this implies an enhancement by a factor 10-30 (depending on the diffusion model adopted) 
of the maximum achievable energy (comparison of the solid, dashed or dotted curves with 
the dot-dashed one). On the other hand, p ma x does not strongly depend on the detailed 
diffusion model adopted, in terms of the turbulence spectrum, with the difference between 
Bohm diffusion and resonant scattering with self-generated Alfven turbulence (comparison 
between the solid and dashed curves) being a factor of 2 at most. A few remarks are worth 
being made. Note that the values of the parameters adopted throughout our calculations are 
reminiscent of those appropriate to the case of young supernova remnants (SNRs). This is 
because we wished to address here the issue of the maximum energy that can be achieved in 
this class of sources, as candidate sources of galactic cosmi c rays. Moreov er there seems to 
be now solid evidence for amplified magnetic fields in SNRs (jWarrenl (120051 )). so that there is 
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Figure 2. We plot the maximum momentum the particles can achieve within a time r = 10 3 yr as a function of the Mach 
number of the shock. The different curves refer to different choices of the diffusion coefficient: the solid line is for self-generated 
diffusion; the dashed line is for Bohm diffusion in the self-amplified magnetic field; the dotted line is for Bohm diffusion in a 
magnetic field that is spatially constant in strength and equal to 8B(0), the value reached at the shock location by the cosmic 
ray streaming induced field; finally the dot-dashed line is for Bohm diffusion in the background magnetic field. 



observational support for the idea that cosmic rays might be efficiently accelerated at SNRs' 
shocks. The shock modification is a natural consequence of the effectiveness of these shocks 
as cosmic ray accelerators. 

The maximum energy that the accelerated particles can achi eve in the context of modified 
shocks is always less than one would naively guess by using a lLagage fc Cesarskyl (jl983al ) 
type of estimate with the upstream velocity set equal to u and the magnetic field evaluated 
at the shock surface. This is due to the combination of two effects. On one side, the net effect 
of the shock modification is to slow down the upstream fluid (precursor), so that the effective 
velocity felt by the particles of given momentum is always lower than uq, thereby making 
the acceleration process slower. This is illustrated effectively through the dotted line in Fig. 
El where the diffusion coefficient is taken as Bohm-like calculated with the magnetic field at 
the shock, assumed constant in the precursor. In this case the value of p max is always lower 
than ~ 10 6 GeV/c for Mach number Mq > 100, while in the absence of shock modification 



12 P. Blasi, E. Amato and D. Caprioli 

one would get p ma x ~ 6 x 10 6 GeV/c. This factor > 5 is solely due to the existence of a 
precursor. There is however a second (less important) effect that causes a further decrease of 
Pmax- the spatial distribution of the accelerated particles in the precursor also determines a 
spatial profile of the strength of the magnetic field associated with the self-generated Alfven 
waves, so that the effective diffusion coefficient is also larger than DBohm{0,p)- This effect 
is visible when comparing the dashed line with the dotted line in Fig. [2j The combination 
of the two effects lead s to an overall decre a se of V rnnx by a factor ~ 10 compared with the 
naive expectation a /a lLagage fc Cesarskyi (1l983al jbl). more evident for higher values of the 
Mach number. 



5 A CRITICAL DISCUSSION OF THE METHOD AND SUMMARY OF 
THE MAIN RESULTS 

In this section we summarize the reasons for the need of the type of calculations presented 
here and we discuss some shortcomings of the proposed mathematical approach. 

To our knowledge this represents the first attempt to determine the maximum momentum 
of the particles accelerated at a cosmic ray modified shock, where the nonlinear dynamical 
reaction of the accelerated particles changes the shocked fluid, the particles self-generate 
their own scattering centers, and the modified system in turn determines the spectral shape 
and energy content of the accelerated particles' distribution. The creation of a shock precur- 
sor in the upstream region and the amplification of the background magnetic field through 
streaming instability are the main reasons to expect that the maximum achievable momen- 
tum may be appreciably different from the one calculated in linear theory for the same value 
of the shock Mach number. The precursor, due to the gradual slowing down of the upstream 
plasma would by itself be expected to slightly decrease the maximum momentum: this effect 
can in fact be seen in Fig. [2} where the dash-dotted line shows the maximum momentum 
as a function of the Mach number when the background field is not amplified. The reduced 
velocity in the upstream plasma slows down the process of return to the shock for the par- 
ticles in the upstream section of the fluid. In the presence of magnetic field amplification 
the maximum momentum is expected to increase. This trend continues up to Mach num- 
bers for which the pressure in the form of accelerated particles cannot increase any more 
(saturation), so that the self-generated magnetic field also saturates. On the other hand the 
compression factor (shock modification) keeps increasing, so that the corresponding p max in 
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Fig. [2] shows a slightly decreasing trend (roughly the same that one obtains in the presence of 
shock modification but without field amplification, dash-dotted line). One should be careful 
not to mistake this trend for a temporal trend. The behaviour of the maximum momentum 
is to be interpreted only as a trend as a function of the fluid Mach number at upstream 
infinity. In a more realistic situation the Mach number can also change with time, which 
might introdu ce additional effects , not ta ken into account here. On the other hand, if the 



conclusions of 



Lagage fc Cesarskyi (If 983bl ) remain valid for this more complex scenario, one 



should expect that most of the acceleration of cosmic rays occurs between the end of the free 
expansion phase and the beginning of the Sedov phase. This is the phase we concentrated 
our attention upon, so that the assumption of temporally constant Mach number adopted 
here is justified. 

The most important equation derived in this work is Eq. [25] for the acceleration time in 



the presence of shock modification. It is important to realiz e that although the ca 



Amato fc Blasil (12005 



culations 



2006|) for 



shown here have been carried out by adopting the theory of : 
the evaluation of the modification and field amplification, Eq. [25] can be applied in a very 
general way. In fact all the quantities involved in the evaluation of Eq. [25] can be determined 
using very different approaches, for instance Monte Carlo simulations, while Eq. [25] can 
be adopted to estimate the instantaneous value of p ma x- Moreover, since Eq. [25] directly 
contains the diffusion coefficients and the dynamical quantities (e.g. fluid velocity), it can 
be applied independently of the specific mechanis m of wave p roduction. If the non-resonant 



wave generation mechanism recently proposed by iBelll (12004 ) were to be responsible for the 
wave production and particle scattering, this would change the mathematical form of the 
diffusion coefficient, and possibly the dynamical reaction of the plasma (if the amplified field 
reaches a dynamically important value), but Eq. [25] could still be applied unchanged. 

The most relevant shortcoming of the calculation illustrated here is in the fact that we 
used the assumption of quasi- st at ionarity in the upstream fluid (not downstream). This 
assumption is fundamental to allow us to apply the method of Laplace transform, which 
can be used only for linear differential equations. At the same time quasi-statio narity is 
also required in order to use our previous calculations (j Amato fc Blasil (120051 . 120061 )) for the 
shock modification and field amplification. These are actually two different aspects of the 
same problem of this approach, which however cannot be avoided if an analytical attempt 
has to be made to calculate the maximum momentum in these complex situations. 

It should however be stressed that the shortcomings just discussed also apply to the 
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approach of lLagage fc Cesarskyl (jl983al ibl): their calculations are only valid for linear shocks 
(Pcr *C Pou%) so that the velocity and density profiles at the shock are spatially constant and 
given. On the other hand the diffusion coefficient, which is self-generated in their approach, is 
time-dependent because the particle distribution function is time-dependent. Despite this, 
the authors fix the diffusion coefficient (and the distribution function of the accelerated 
particles) and proceed by applying the method of the Laplac e transform. In this sense o ur 



Lagage &: Cesarskyl (ll983al Jbh 



approach should be considered as the extension of the work of 
to the case in which the dynamical reaction of the accelerated particles is not neglected. 
The common finding of several independent approaches to particle acceleration at shocks 
is that shock modification manifests itself in a prominent way, as a consequence of efficient 
acceleration, therefore our extension seems particularly important at this time. 

We conclude this section by providing a short, itemized summary of our main scientific 
results: 

1) the maximum momentum of accelerated particles can be as hi gh as ~ 10 6 GeV/c 



for pro tons, when all effects are taken into account. The approach of lLagage fc Cesarsky 
( 1l983al jbl) would give p ma x ~ 6 x 10 6 GeV/c when the magnetic field is assumed spatially 
constant to its maximum value at the shock surface. The same approach for 6B /B w 1 would 
reproduce the well known value p ma x = 10 5 GeV/c for B Q = IO/iC and u = 5000km/s. 

2) the net effect of the shock modification is to create a precursor in both the fluid 
velocity and magnetic field strength. The former decreases while moving towards the shock 
while the latter increases. Both phenomena contribute to lowering p max with respect to the 
value ~ 6 x 10 6 GeV/c. 

3) Provided the magnetic field is amplified by the streaming instability of cosmic rays, 
the value of p ma x is weakly dependent upon the exact spectrum of the magnetic turbulence. 

4) The expression provided in Eq. [25] is the only expression presently available that allows 
one to calculate the acceleration time for cosmic ray modified shock waves. This expression 
is very general and can be applied independently of the specific treatment of the shock 
modification and reduces to the well known formula for weakly modified shocks. 
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